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ABSTRACT 

Recent spectropolarimetric studies of 7 SPB and (3 Cep stars have suggested that photospheric 
magnetic fields are more common in B-type pulsators than in the general population of B stars, 
suggesting a significant connection between magnetic and pulsational phenomena. We present an 
analysis of new and previously published spectropolarimetric observations of these stars. New Stokes 
V observations obtained with the high-resolution ESPaDOnS and Narval instruments confirm the 
presence of a magnetic field in one of the stars (e Lup), but find no evidence of magnetism in 5 
others. A re-analysis of the published longitudinal field measurements obtained with the low-resolution 
FORSl/2 spectropolarimeters finds that the measurements of all stars show more scatter from zero 
than can be attributed to Gaussian noise, suggesting the presence of a signal and/or systematic under- 
estimation of error bars. Re-reduction and re-measurement of the FORSl/2 spectra from the ESO 
archive demonstrates that small changes in reduction procedure lead to substantial changes in the 
inferred longitudinal field, and substantially reduces the number of field detections at the 3a level. 
Furthermore, we find that the published periods are not unique solutions to the time series of either 
the original or the revised FORSl/2 data. We conclude that the reported field detections, proposed 
periods and field geometry models for a Pyx, 15 CMa, 33 Eri and V1449 Aql are artefacts of the 
data analysis and reduction procedures, and that magnetic fields at the reported strength are no more 
common in SPB//3 Cep stars than in the general population of B stars. 

Subject headings: stars: magnetic — stars: pulsators — stars: individual(a Pyx, 15 CMa, 33 Eri, e 
Lupi, HY Vel, V1449 Aql) 



1. INTRODUCTION 

Magnetism in hot, massive stars appears to be a rela- 
tively uncommon phenomenon. Despite tremendous ob- 
servational resources directed to large-scale surveys since 
the beginning of the 21st century only a few dozen mag- 
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netic OB stars have been firmly identified. Their rarity 
notwithstanding, it is apparent that those stars which 
do host detectable fields have certain common magnetic 
characteristics: the fields are usually topologically dipo- 
lar, with typical inferred strengths of hundreds to thou- 
sands of gauss. There often exists a significant obliquity 
between the magnetic and rotational axes. The mea- 
sured fields are remarkably stable, persisting over many 
rotational cycles with no detectable secular change of 
field strength or geometry. The stars hosting these fields 
could not, however, be more diverse: they are both old 
and young, with strong and weak winds, and rotational 
periods varying from less than a day to decades. Some 
possess photospheric chemical peculiarities, others winds 
or circumstellar matter that interacts with the field. 0th- 
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crs experience pulsations. 

While magnetic OB stars generally provide fertile 
ground for exploring important facets of stellar evolu- 
tion, the pulsating stars are of particular interest. Pul- 
sating stars provide a unique opportunity to study the 
subtle influence of the magnetic field on changes of stel- 
lar structure in near real-time. Pulsating stars also offer 
the potential to probe their interiors via asteroseismol- 
ogy, yielding information about stellar structure. Pul- 
sating stars have furthermore been at the forefront of 
recent discussions of the incidence of magnetism in hot 
stars. Hubrig et al. (2006, 2009) reported magnetic 
measurements of 61 pulsating or candidate pulsating B- 
type stars (50 slowly pulsating B (SPB) stars aud 19 /3 
Cep stars), acquired using the FORSl spectropolarime- 
ter at the ESO VLT. In this large sample, Hubrig et al. 
(2009) claimed 3a detection of 52 percent of the SPB 
stars (i.e. 26 of 50 SBP stars) and 31 percent of the /3 
Cep stars (6 of 19 stars). These results suggested that 
magnetic fields are significantly more common in pulsat- 
ing B-type stars than in the overall population of B-type 
stars. However, Silvester et al. (2009), analysing inde- 
pendent magnetic observations of pulsating B-type stars 
obtained with the ESPaDOnS and Narval spectropo- 
larimeters, found a much lower incidence of magnetic 
stars. In particular, of 8 stars claimed to be magnetic 
by Hubrig et al. (2006, 2009) that were common to both 
the Hubrig et al. and Silvester et al. studies, Silvester et 
al. (2009) confirmed fields in only 2. Furthermore, the 
Magnetism in Massive Stars (MiMeS) survey component 
has found magnetic fields to be no more common in SPB 
or (3 Cep stars than in B stars in general (Grunhut et al., 
2011). These represent serious discrepancies between the 
results of these two different approaches (investigated in 
detail by Bagnulo et al., 2012). 

Recently, Hubrig et al. (2011a; hereafter H2011a) re- 
ported numerous longitudinal magnetic field measure- 
ments of 6 SPB and ,5 Cep stars, obtained from low- 
resolution (A/AA ~2000) spectropolarimetric observa- 
tions acquired using the FORSl and F0RS2 instruments 
at ESO's Very Large Telescope (VLT). All six program 
stars - 33 Eri (HD 24587, B5V, SPB), CMa (HD 
46328, B0.7IV, /3 Cep), 15 CMa (HD 50707, BlIV, 
13 Cep), a Pyx (HD 74575, B1.5IH, /3 Cep), HY Vel 
(HD 74560, B3IV, SPB), and e Lupi (HD 136504, B2IV, 
$ Cep) - were reported by those authors to host photo- 
spheric magnetic fields with maximum measured longitu- 
dinal components of order 100 G. Models of the siirface 
magnetic field geometry in the context of the Oblique 
Rotator Model (ORM; including rotational periods Prot, 
rotational axis inclination angles i, magnetic obliquity 
angles /3, and dipolar field strengths B^) were reported 
for 4 of the 6 stars - 33 Eri, CMa, 15 CMa, and a 
Pyx. The study by H2011a was followed by a similar 
examination by Hubrig et al. (2011b, hereafter H2011b) 
of another /3 Cep star, V1449 Aql, based on observa- 
tions obtained with the SOFIN spectropolarimeter on 
the Nordic Optical Telescope. H2011b report an appar- 
ently coherent sinusoidal modulation of the longitiidinal 
field with a full amplitude of ^1 kC. These results are 
qualitatively distinguished from most previous studies of 
magnetism in pulsating stars by the long-term monitor- 
ing of a large sample of particular stars. They potentially 



represent an important contribution to the accumulated 
knowledge and scientific discussion regarding magnetism 
of pulsating massive stars. 

The aim of the current paper is to further explore the 
incidence of magnetic fields in pulsating B-type stars by 
independently investigating and elaborating the results 
reported by H2011a and H2011b. In Section 2 we sum- 
marise the observational data employed in this study: 
new ESPaDOnS and Narval observations, the published 
FORSl/2 longitudinal magnetic field measurements of 
H2011a, and new magnetic field measurements obtained 
from re- reduction and re- measurement of the FORSl/2 
spectra acquired by H2011a (since the SOFIN observa- 
tions reported by H2011b remain proprietary, we are 
unable to perform as comprehensive an analysis in this 
case). In Section 3 we present the magnetic analysis of 
the ESPaDOnS and Narval spectropolarimetry, the orig- 
inal published FORSl/2 results, and the measurements 
from re-reduced FORSl/2 observations extracted from 
the ESO archive. In Section 4 we compare the results 
between instruments and reductions. In Section 5 we 
evaluate the models presented by H2011a and H2011b. 
In Section 6 we discuss the results and their implica- 
tions for the magnetic fields of the individual targets, the 
general incidence of magnetic fields in pulsating B-type 
stars, and the performance and reliability of the FORSl 
and F0RS2 instruments as stellar magnetometers. 

2. OBSERVATIONS 
2.1. ESPaDOnS and Narval 

High-resolution (R^65000) circular polarization 
(Stokes V) spectra were obtained for 6 of the stars 
studied by H2011a and H2011b using the cross-dispersed 
echelle spectropolarimeters ESPaDOnS (at the Canada- 
France-Hawaii Telescope) and its clone Narval (at 
the Telescope Bernardot Lyot (TBL), Pic du Midi 
Observatory) within the context of the MiMeS Large 
Program. Altogether, 18 independent observations were 
obtained between February 2010 and June 2011. Each 
spectropolarimetric sequence consisted of four individ- 
ual subexposures taken in different orientations of the 
half-wave retarders. From each set of four subexposures 
we derived Stokes / and Stokes V spectra following the 
double-ratio procedure (e.g. Donati et al., 1997). All 
frames were processed using the automated reduction 
package Libre-ESpRIT (Donati et al., 1997), with 
the Stokes / spectra normalized order-by-order using 
Legendre polynomials following extraction of Stokes V. 
The peak signal-to- noise ratios (SNRs) per 1.8 km s^^ 
spectral pixel in the reduced spectra range from about 
450 to about 1300. 

As discussed by Semel et al. (1993) and Schnerr et al. 
(2006), rapid temporal changes in line profile shapes and 
positions due to pulsation may also introduce spurious 
signal in the inferred polarisation spectra. In the case 
of a non-magnetic star, such a phenomenon may result 
in a false magnetic detection. Spurious signatures may 
equally be introduced by pulsation when a magnetic field 
is present, combining with the real Zeeman signature and 
potentially modifying its shape and amplitude. 

Since all targets are either confirmed or candidate /3 
Cep or SPB stars, subexposures were timed, as far as 
possible, to correspond to no more than ~l/80th of a 
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pulsation cycle to ensure that no significant variability 
occurred during the observation. Of the five observed 
stars for which periods are known, three (e Lupi (pul- 
sation periods of 0.155 d (primary) and 0.096 d (sec- 
ondary); Uytterhoeven et al. 2005, Telting et al. 2006), 
33 Eri (0.864 d; de Cat et al, 2005) and HY Vel (3.102 d; 
Lefevre et al., 2009)) were observed with subexposures 
of hundredths of a period. 15 CMa (0.093 d - 0.181 
d; Shobbrook et al., 2006) was observed with subexpo- 
sures corresponding to less than 1 / 80th of the longer pe- 
riod, but about l/60th of the shorter period. V1449 Aql 
is a multiperiodic star with a dominant mode period of 
0.182 d (Briquet et al., 2009) and several lower-amplitude 
modes as rapid as 0.119 d. Longer subexposures were 
required due to its faint apparent magnitude, resulting 
in subexposures for this star that are somewhat longer 
than l/80th of these modes, a Pyx is a candidate /? Cep 
star for which the period is not yet known. However, its 
subexposure durations were 150 s, and therefore would 
satisfy the limit described above as long as its pulsation 
period is longer than 0.14 d. 

Given the uncertainties described above, we have per- 
formed a number of tests to establish that pulsations did 
not significantly impact the magnetic diagnosis. We ex- 
amined the Stokes / profiles associated with each of the 
subexposures of each observation. For e Lupi, a Pyx, 
33 Eri, 15 CMa and HY Vel, we confirm that no vari- 
ation of the / profile during an observation is detected 
above the noise. For V1449 Aql, on the other hand, the 
/ profile is observed to vary during an observation by a 
maximum of about 10 km s~^ in radial velocity, with a 
small (~ 10%) variation in depth, and negligible change 
in shape. To verify that instrumental instabilities or line 
profile variations such as those observed in V1449 Aql do 
not translate into a spurious Stokes V signal, we com- 
puted diagnostic null polarization spectra (labeled K) 
by combining the four subexposures in such a way that 
stellar polarization cancels out (see Donati et al. 1997 
or Bagnulo et al. 2009 for more details on the definition 
of . In no case is any signal detected in N above the 
level of the noise. Considering the numerical simulations 
of Schnerr et al. (2006) which demonstrate that spurious 
signatures due to pulsation appear in both V and N , this 
confirms in the case of V1449 Aql that the observed line 
variations introduce no spurious signal into Stokes V . 

Notwithstanding this conclusion, variability of the pro- 
file of V1449 Aql during the observation may result in 
temporal smearing of the Stokes / and V profiles. To 
verify that this did not significantly reduce our sensitiv- 
ity to any magnetic field, we have performed a test in 
which we generated fictitious non-zero Stokes V profiles 
by computing the numerical derivative of the Stokes / 
profile corresponding to each of the subexposures. While 
these profiles do not correspond to any real magnetic 
field (e.g. any prediction of the model of H2011b), they 
have the qualitative characteristics of a real Stokes V 
Zeeman signature, and in particular reflect the changing 
shape and position of the line profile of V1449 Aql during 
the course of a polarimetric observation. They therefore 
serve to understand the impact of the temporal smear- 
ing on Stokes V . The fictitious V profiles exhibit small 
changes from subexposure to subexposure, refiecting the 
radial velocity, depth and shape variability of Stokes /. 
Nevertheless, we find that the characteristics of the mean 



V profile inferred by averaging the subexposures does 
not differ qualitatively from the V profiles of the indi- 
vidual subexposures. In particular its overall shape and 
amplitude are not significantly modified. Therefore tem- 
poral smearing of the ESPaDOnS and Narval Stokes V 
and / profiles should not significantly affect either the 
detectability or the interpretability of any magnetic sig- 
natures present in the spectrum of V1449 Aql. 

A related secondary effect involves the depth- 
dependence of the pulsational dynamics, and therefore 
a (weak) potential systematic difference in the shape of 
the line profiles as a function of line strength. To evalu- 
ate whether the depth of formation of spectral lines has 
any infiuence on the magnetic diagnosis of V1449 Aql, 
we performed a simple test using Least Squares Decon- 
volution (LSD; described in detail in Section 3.1) masks 
filtered for weak (unbroadened central depth < 0.2 of the 
continuum) and strong (> 0.2 of the continuum) spectral 
lines. The two masks contain similar numbers of spectral 
lines. We used these masks to extract Stokes V profiles 
of e Lupi, the only star of the sample in which a magnetic 
field is detected. The LSD Stokes / profiles extracted us- 
ing the two masks are reasonably similar in shape; the 
asymmetric structure of the line profile (see Fig. 1) is 
marginally clearer in the "weak" line, while the wings of 
the "strong" line are more extended than those of the 
"weak" line. The two LSD profiles differ by about a fac- 
tor of 1.5 in depth. Longitudinal fields measured from 
the "strong" and "weak" line LSD profiles both corre- 
spond to magnetic detections, and agree within the error 
bars. 

While the two tests described above do not represent 
a sophisticated evaluation of the temporal smearing and 
depth-dependence effects on magnetic diagnosis, they il- 
lustrate that no significant effect is present at the preci- 
sion of the ESPaDOnS and Narval observations described 
in this paper. 

In this paper, we present 4 observations of 33 Eri, 4 
of 15 CMa, 3 of a Pyx, 2 of e Lupi, 4 of V1449 Aql 
and 1 of HY Vel. In the case of HY Vel, poor observing 
conditions limited the quality of the data and prevented 
us from coming to any firm conclusions. 

The log of CFHT and TBL observations is reported 
in Table [T] For e Lupi and V1449 Aql, co-added spec- 
tra were created from observations taken on the same 
night in order to increase the SNR. In the case of V1449 
Aql radial velocity differences (described above) between 
observations were significant enough that spectra were 
brought to a common velocity before being combined. 

2.2. FORSl/2 and SOFIN 

We have re- analyzed the FORSl/2 observationtQ ob- 
tained by H2011a, both in published and re- reduced 
forms. H2011a presented 62 Stokes V observations of 5 of 
our program stars, obtained between November 2003 and 
March 2010 using grism 600B, delivering a spectral res- 
olution R~2000. The observations sample this time pe- 
riod nonuniformly, and the number of observations varies 
significantly from star to star (7 observations each for 
a Pyx and e Lup, 15 observations each for 15 CMa and 
HY Vel, and 18 observations for 33 Eri). 

^ After decommissioning in March 2009, the polarimetric optics 
of FORSl were moved to FORS2. 
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TABLE 1 

Summary of results from ESPaDOnS and Narval data. Note that subexposure times correspond to 1/4 of the total 

exposure time. 
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We also analyzed the 2 FORSl/2 and 13 SOFIN mea- 
surements of V1449 Aql reported by H2011b. The 
SOFIN spectra were obtained from August-September 
2009 and in July 2010, and were made at intermediate 
spectral resolution (R^30,000). 

3. MEASUREMENTS OF LONGITUDINAL MAGNETIC 
FIELDS 

3.1. Measurements from ESPaDOnS and Narval spectra 

The ESPaDOnS/Narval spectra were analyzed using 
the Least Squares-Deconvolution (LSD) multiline analy- 
sis method developed by Donati et al. (1997). LSD com- 
bines the information from essentially all metallic and He 
lines in the spectrum by means of the assumption that 
the spectrum can be reproduced by the convolution of a 
"mean" line profile (the "LSD profile" ) with an underly- 
ing spectrum of unbroadened atomic lines of given line 
depth, Landc factor, and wavelength (the "line mask" , 
as described by Wade et al., 2000). This process allows 
the computation of averaged Stokes /, V and N profiles 
with much higher signal-to-noise ratios than those of the 
individual lines, dramatically improving the detectabil- 
ity of Zeeman signatures due to stellar magnetic fields. 
The application of LSD to spectra of pulsating B stars is 
discussed in detail by Silvester et al. (2009). 

The line masks used in this analysis were modified from 
generic masks obtained using the Vienna Atomic Line 
Database (VALD; Piskunov et al., 1995; Ryabchikova et 
al., 1997; Kupka ct al., 1999; Kupka ct al., 2000) ex- 
tract STELLAR requests. The VALD line lists arc char- 
acterized by the stellar effective temperatures and grav- 
ities of the individual stars (as reported by H2011a), so- 
lar abundances, and a minimum line depth of 10% of the 
continuum. The masks were customized to the spectrum 
of each star by interactive comparison, removing weak 
linos and those blended with Balmcr or telluric lines. 
This 'cleaning' of the mask ensures that lines included 
in the LSD procedure have shapes that are as similar 
and uniform as possible. 

The lines remaining in the mask following the cleaning 



procedure - between 230 and 320 lines depending on the 
physical and spectral properties of the star - were then 
adjusted to better match the observed line depths in the 
associated observed spectnun. 'Tweaking' the lines in 
this fashion improves the fit between the actual spectrum 
and the LSD spectrum model (that is, the convolution 
of the LSD profile and the line mask). 

Two methods were used to evaluate the presence of a 
photospheric magnetic field: a statistical test performed 
on the Stokes V profile, and direct inference based on 
the significant detection of a longitudinal magnetic field. 
The former method, described by Donati, Semel, &: Rees 
(1992) and Donati et al. (1997), employs the reduced 
of the signal in Stokes V within the botmds of the 
line profile. It reports the detection of a magnetic sig- 
nature as 'definite' if the formal detection probability is 
greater than 99.999%, while the detection probabilities 
outside the Stokes V line profile, and inside the diagnos- 
tic null {N) line profile, are both negligible; a 'marginal' 
detection corresponds to a detection probability between 
99.9% and 99.999%. 

The second, direct inference method involves compu- 
tation of the mean longitudinal magnetic field from the 
first-order moment of the Stokes V profile within the line 
according to the expression: 



-2.14 X 10^ 



/ vV{v)dv 



>^9esc J [7c - I{v)] dv 



(1) 



where (Bz) is the mean longitudinal magnetic field in 
G, V is the velocity in km s^^ within the profile mea- 
sured relative to the centre of gravity (Mathys et al., 
1989; Donati et al. 1997; Wade et al., 2000), and A 
and ^eflF Siice the reference values of the wavelength in nm 
and Lande factor used in computing the LSD profiles. 
This equation is also applied to the null profile, yield- 
ing a comparable measurement (Nz). The uncertainties 
associated with (Bz) and (Nz) wore determined by prop- 
agating the formal (photon statistical) uncertainties of 
each pixel through Eq. (1). This is described in more 
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Fig. 1. — Examples of LSD profiles of stars observed with ESPaDOnS and Narval. Stokes / in black (bottom), diagnostic N in blue 
(middle), V in red (top). Both V and A' have been scaled and shifted for display purposes. Vertical lines denote the integration range 
adopted for longitudinal field measurement. Top, left-right: e Lupi, V1449 Aql, a Pyx, Bottom, left-right: 33 Eri, 15 CMa, HY Vel. Note 
the similarity between A'' and in 5 of the 6 stars. Note also that the scaling of V and A'' varies substantially between stars. For each 
star, the LSD profile with the highest SNR is shown. 



detail by Silvester et al. (2009). 

The integration ranges employed in the evaluation of 
Eq. (1) associated with each LSD profile were selected 
individually through visual inspection so as to include 
the entire span of the Stokes / profile. These ranges are 
indicated in Fig. [TJ and reported for each observation 
in Table [TJ Integration ranges can differ significantly for 
different observations of the same star, especially when 
pulsations are a factor as they are here, and the ability to 
tailor the integration range to the shape of each line pro- 
file and thus optimize the magnetic diagnosis is an im- 
portant strength of high-resolution spectropolarimetry. 
The pulsational character of the stars is clearly visible in 
the often rather asymmetric LSD Stokes / profiles. We 
tested the sensitivity of the inferred longitudinal field 
and its error bar to the integration range by adjusting 
the integration limits by ±10%. This has the effect of 
decreasing (for smaller range) or increasing (for larger 
range) the error bar by 15-20%, but otherwise has no 
impact on the results. These results are qualitatively 
consistent with experiments carried out by Neiner et al. 
(2012) (see their Fig. 3). 

Most stars show relatively broad wings, attributable 
to the inclusion of He lines in the line masks. Remov- 
ing the He lines from the mask reduces this effect at the 
expense of sensitivity, yielding an increase in the error 
bars of {Bz) on order of 50%. Experiments using e Lupi 
(the only star of the sample in which magnetic field is 
detected in the ESPaDOnS/Narval spectra) clearly show 
that masks including and excluding He lines yield congru- 
ent results in the presence of a magnetic field (apart from 
differences in error bars). The line masks from which the 
LSD profiles shown in Fig. [T] were constructed include 
He lines. 

3.2. Published measurements 



The mean longitudinal magnetic field (Bz) of FORSl/2 
data was derived using the linear regression method de- 
veloped by Bagnulo et al. (2002). The measurements 
were obtained in two ways: using only hydrogen Balmer 
lines and using the full spectrum (including H Balmer 
lines along with metal and He lines). No diagnostic null 
measurements were reported. We assume that uncertain- 
ties were derived from the formal uncertainties obtained 
from the regression. These measurements, which we re- 
analyze below, are reported in Tables 2 and 3 of H2011a. 

(Bz) was measured from SOFIN using the moment 
technique of Mathys (1994). Again no diagnostic nulls 
are reported. The measurements are provided in Table 
1 of H2011b. 



3.3. Measurements from re-reduced spectra 

In order to test the robustness of the field measure- 
ments reported by H2011a, the archival data from which 
they were determined were obtained and re-reduced us- 
ing a new FORS reduction pipeline, the details of which 
are described by Bagnulo et al. (2012). We have mea- 
sured the longitudinal field using the method of Bagnulo 
et al. (2002) from H lines (H/3 to the Balmer jump), 
from the full spectrum (from 3800 A-4950 A) , as well as 
from "metal lines" (i.e. the full spectrum excluding H 
lines) . The measurements were performed using a larger 
aperture for spectrum extraction, and omitting the re- 
gion 3900-4000 A, which has been found to be afflicted 
by internal refiections. For all three sets of measurements 
obtained from the Stokes V spectra, longitudinal fields 
from the corresponding diagnostic null spectra were also 
measured. Full spectrum and H line measurements are 
compared to the original published measurements in Fig. 
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a Pyx 





Fig. 2. — Longitudinal magnetic field measurements obtained with FORSl/2 (open blue triangles, from H2011a and H2011b) and 
ESPaDOnS/Narval (filled red diamonds). The data have been phased according to the ephemerides proposed by H2011a and H2011b. All 
ESPaDOnS/Narval measurements are consistent with the absence of magnetic field: the dotted line shows (Bz) = 0, while the solid line 
shows the field curve from the H2011a and H2011b models. Top-Bottom: a Pyx, 33 Eri, 15 CMa, V1449 Aql 

in the data: 12 full spectrum measurements and 10 H 
line measurements are non-zero at greater than Scr sig- 
nificance. Moreover, comparison of the measurements 
for each star with a variety of simple models consistently 
favours the presence of a variable magnetic field. On 
the other hand, the measurements from the re-reduced 
FORSl/2 spectra present far fewer significant detections, 
and remarkably these correlate rather poorly with the de- 
tections in the original analysis. Clearly the characteris- 
tics of the new measurements differ from those previously 
published. 

The published FORSl/2 measurements for all stars are 
compared to the re-reduced values in Fig. [3l Visual 
inspection of Fig. [3] shows that while differences exist 
between individual (Bz) measurements, the longitudinal 
fields determined from the two data sets are generally 
in formal agreement. However, the error bars of the 
re-reduced data are in many cases systematically larger 
than those of the original published results, leading to 
a greater proportion of the individual measurements be- 
ing formally consistent with a null result; this may be 
because the error bars of the re-reduced measurements 
are weighted with the square root of the reduced x^- 
Nevertheless, the re-reduced measurements still present 



proprietary and so a re-reduction cannot be undertaken. 
4. RESULTS 

The results of the LSD analysis are summarized in Ta- 
ble [U and illustrated in Figs. [T] and [2] According to 
our detection criteria, just one of the 6 stars - e Lupi - 
shows evidence for the presence of a magnetic field in its 
photosphere, e Lup exhibits a significant Stokes V LSD 
profile, corresponding to longitudinal magnetic fields de- 
tected at 5.3(7 and 4. Scr. The remaining five stars exhibit 
no significant longitudinal magnetic field nor Stokes V 
signal. 

To evaluate the robustness of the magnetic field diag- 
nosis from Stokes V, we have also computed both the 
detection probability and longitudinal field using the N 
LSD profiles. In no case is any significant polarised signal 
or longitudinal field detected in the those data. 

The remarkable result of the ESPaDOnS and Narval 
observations is the absence of any evidence of magnetic 
field for four stars - a Pyx, 15 CMa, 33 Eri and V1U9 
Aql - for which magnetic field geometry models were pre- 
sented by H2011a and H2011b. 

Examination of the published measurements of a Pyx, 
HY Vel, 33 Eri and 15 CMa suggest that signal is present 
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Fig. 3. — Comparison between FORSl/2 data reductions. Left: 
full spectrum measurements; right: H Balmer line measurements. 
The straight line corresponds to x = y. A magnetic field is unam- 
biguously detected in both reductions in e Lupi only; in all other 
cases the re-reduced data are most consistent with the originally 
published data around {B^) = 0. In the online version, 3cr out- 
liers are indicated in color: green for 3(t measurements reported 
by H2011a, blue for Scr outliers in the re-reduced measurements, 
and red for measurements significant in both reductions. Note that 
some such outliers agree in significance but not polarity. 



a number of 3a detections that are difficult to reconcile 
with Gaussian noise in the absence of a magnetic field. 



It is hard to overstate the importance of examining 
diagnostic nulls as well as circular polarization spectra 
when examining spectropolarimetric data. In the ideal 
case, there would be either a clear signature in Stokes V 
and no signature in N (an unambiguous detection), or 
no signature in Stokes V and no signature in A'^ (an un- 
ambiguous non-detection.) However, it occasionally hap- 
pens that there is a signature in both Stokes V and N, 
or in but not in Stokes V. Field detections in the null 
profiles may be due to either instrument fiexures and/or 
radial velocity shifts intrinsic to the source. Schnerr et 
al. (2006) have modelled the effects of stellar pulsations 
on both Stokes V and in the null profiles. Bagnulo et al. 
(2012) have argued that magnetic field measurements ob- 
tained with FORS represent a slightly more complex sit- 
uation, because radial velocity shifts affect not only null 
profiles (and, to a lesser extent, also Stokes V profiles), 
but also the normalised derivative of Stokes /. Bagnulo 
et al. (2012) have concluded that while a field detection 
in the null profiles does not necessarily imply that a field 
detection in Stokes V (obtained from the same dataset) 
is spurious, it demonstrates that photon noise is not the 
only contribution to the error bars. Examining the inci- 
dence of longitudinal field detections from the diagnostic 
null spectra of pulsating B stars observed with FORS 
confirms that the error bars published by H2011a (and 
to a lesser degree those obtained in our re-reduction) are 
an underestimate of the actual error bars by a factor of 
30 to 50%. 

While examination of N is indispensable, it is not 
always sufficient. For instance, in the case of a Pyx, 
the FORS spectrum taken on HJD 2454109.150 yields a 
strong (Bz) measurement in the blue end of the spec- 
trum but shows nothing in the Balmer lines II/3-H5 or in 
the metal lines. When (Bz) measurements are extracted 
separately for H/3-H(5 and the higher Balmer lines, the 
former group yields (B^) ~ —80 ± 150 G from Stokes V 
and ~ 130 ± 110 G for N (consistent with a null result), 
while the latter set of lines yields (B^) ~ -1900 ±260 G 
from Stokes V and ~ —625 ±200 G from N, an apparent 
7.3(7 detection (with a 3.1ct detection in A^). This seems 
unlikely to happen if the circular polarization signature 
arises from real photospheric magnetic fields. 

5. EVALUATION OF MODELS 
5.1. ESPaDOnS and Narval 

Recall that of the 7 pulsating B stars claimed to be 
magnetic by 112011a and H2011b, ephemerides are pre- 
sented for 5: CMa, 33 Eri, 15 CMa, a Pyx and V1449 
Aql. For each of these stars they derived a model of the 
magnetic field, including inclination of the rotation axis 
relative to the line of sight i, magnetic obliquity /3 and 
dipolar field strength B^. 

ESPaDOnS and Narval longitudinal magnetic fields of 
15 CMa, 33 Eri, a Pyx and V1449 Aql are uniformly con- 
sistent with the absence of any magnetic fields, in con- 
trast to the FORSl/2 resuhs published by 112011a and 
112011b. One possible explanation of this discrepancy is 
suggested by an examination of Fig. [2] While the ES- 
PaDOnS/Narval {Bz) measurements are all consistent 
with a null longitudinal field, many of these measure- 
ments were obtained at times when, according to the pro- 
posed ephemerides and field geometry models, the longi- 



8 



Shultz et al. 



a Pyx, ^ = Q.^Q2, model 33 Eri, = 0.787, model 15 CMd. = 0.677, mode V1449 Aql, = 0.^29, mode 
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Fig. 4. — Comparison of observed with synthethic LSD profiles. Synthetic LSD profiles above, corresponding observed profiles below; 
left-right: a Pyx, 33 Eri, 15 CMa, V1449 Aql. The models are created according to the parameters of H2011a and H2011b and phased 
according to the ephemerides therein. Turbulence has been added to improve the agreement in the wings of the intensity profile. Synthetic 
Gaussian noise corresponding to that of the observed profiles has been added. In all cases the Stokes V signatures predicted by the proposed 
field configurations should have been clearly detectable at the observed phases and SNRs. 



tudinal field should be close to zero. Some data points, 
because of the relatively low amplitudes of the model lon- 
gitudinal field variations and generally larger error bars 
of the FORS data, are in fact in reasonable agreement 
with the FORS measurements obtained at those same 
phases. So, while some of the ESPaDOnS/Narval (Bz) 
measurements are strongly inconsistent with the pub- 
lished variations - V1449 Aql in particular - one can 
reasonably accommodate many of the remaining differ- 
ences within the measurement uncertainties and possibly 
the uncertainties of the model parameters. 

However, the longitudinal field is but one (secondary) 
diagnostic of the magnetic field extracted from the ES- 
PaDOnS and Narval data. This is because a large vari- 
ety of magnetic configurations can produce a mean lon- 
gitudinal field component that is formally null (i.e. for 
which the first-order moment of Stokes V is zero) at some 
phases. Nearly all of these configurations will generate 
a detectable Stokes V signature in the velocity-resolved 
line profile. The morphology of the Stokes V profile of- 
ten allows these various magnetic configurations - all of 
which yield the same (zero) longitudinal field ~ to be 
distinguished. 

None of the LSD Stokes V profiles for a Pyx, 15 CMa 
33 Eri or V1449 Aql show any significant signal. To pre- 
dict the Stokes V profiles we would expect to observe as 
a consequence of the proposed field geometry models, we 
calculate synthetic disc-integrated Stokes / and V pro- 
files. The model is essentially identical to that described 
by Petit et al. (2008, 2011): the local Stokes / profiles 
are assumed to be Gaussian in shape; the local Stokes 
V profile is computed assuming the weak-field approxi- 
mation; and an explicit disc integration is performed as- 
suming a simple limb-darkening law. Various calculated 
Stokes / and V profiles computed by the simple model 
have been compared to calculations obtained using the 
Zeeman code (Landstreet et al. 1988, Wade et al. 2001) 
and are found to be in acceptable agreement with the 
detailed polarized radiative transfer calculations. 

We proceed to compute "synthetic observations" as fol- 
lows, treating the LSD profiles as a single, mean spectral 
line. First, we match as closely as possible the observed 



and model Stokes / profile. We begin with the vsini re- 
ported by H2011a and H2011b, scale the model profile to 
match the depth of the observed LSD profile, then add 
Gaussian turbulence as necessary to reproduce as closely 
as possible the line wings. Then, using the mean wave- 
length and Landc factor of the real LSD profile, we com- 
pute the model Stokes V profile corresponding to the ap- 
propriate phase and field geometry parameters {i, 13, B^) 
of H2011a and H2011b. Finally, we resample the model 
profiles to the 1.8 km s~^ pixels of the observed profiles 
and add synthetic Gaussian noise corresponding to the 
S/N of the real LSD profile. Clearly, these fits to the LSD 
profiles are not meant to represent a detailed line syn- 
thesis in the complex pulsating atmosphere. However, 
given the success of Donati et al. (2001) in applying 
a similar modeling procedure to self-consistently repro- 
duce the Stokes V profiles and longitudinal field curve of 
/3 Cep itself (which we point out appears to exhibit much 
larger pulsational distortion of its line profiles than the 
stars discussed here) we are confident that this proce- 
dure provides a reasonable first-order calculation of the 
amplitude and shape of the Stokes V profiles predicted 
by the H2011a and H2011b models. 

Analysing the synthetic observations using the same 
procedures as used for the real observations, we find lon- 
gitudinal fields that are (naturally) in good agreement 
with the model predictions (i.e. the solid curves in Fig. 
[2]). The Stokes V profiles of the synthetic observations, 
on the other hand, are strikingly different from those 
observed. In every observation of 33 Eri, 15 CMa and 
a Pyx, the proposed field geometry models predict Zee- 
man signatures that would be easily observed and defi- 
nitely detected i.e. detection probability > 99.999%) in 
our observations. These results are illustrated in Fig. 2] 
for LSD profiles obtained from masks including He lines; 
if He lines are excluded, the results are essentially the 
same. The observations included in this figure were se- 
lected to be those for which the longitudinal field was 
most compatible with the models, as determined from 
Fig. [21 The disagreement between the predicted and 
observed Stokes V profiles is quite striking. Given the 
SNRs obtained in ESPaDOnS and Narval spectra, the 
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Period (days) 

Fig. 5. — Top: Periodograms generated using the originally 
published and re-reduced FORSl/2 data. H2011 period is indi- 
cated with a red arrow. Systematically larger error bars in the re- 
reduction lead to a smaller reduced x^- The periodograms share 
a similar variability pattern, presumably beats arising from the 
window function. The most significant peaks are not at the same 
periods. Bottom: CLEANed periodograms. Once again, the de- 
tected periods are dissimilar between reductions. 

models proposed by H2011a and H2011b predict definite 
detections of Zeeman signatures in all 15 observations of 
33 Eri, 15 CMa, a Pyx and V1449 Aql. Therefore, unlike 
the longitudinal field measurements, the incompatibility 
of the observed Stokes V profiles and those predicted by 
the models is totally unambiguous. 

5.2. FORSl/2 and SOFIN 

One puzzling aspect of the results presented by H2011a 
and H2011b is the remarkably low values of the reduced 
for some stars when the data are fit with a sinusoidal 
variation: in general a reduced significantly below 1 
for a sample of this size is indicative of fitting noise or 
over-estimation of error bars. However, we recall that 
the re-reduction of the FORSl/2 spectra has yielded er- 
ror bars systematically larger than those of the original 
reduction, suggesting that the H2011a error bars are not 
substantially over-estimated (and are potentially under- 
estimated.) 

To investigate this matter further, we have constructed 
periodograms using both the full spectrum data and the 
H line data, both published and re-reduced, for all five 
stars for which detailed models were presented by H2011a 
and H2011b. We use a Lomb-Scargle type approach, 
computing the reduced as a function of adopted pe- 
riod by fitting a linear first-order sinusoid to the data. 
Examples of the derived periodograms are shown in Fig. 
[5l The periodograms are quite complex, and in all cases 
numerous periods exist corresponding to fits to the data 
with reduced x^s comparable to the periods reported by 
H2011a and H2011b. In other words, the H2011a and 
H2011b periods are not unique. None of these periods 
seems to be a priori preferable to another. In the cases 
of a Pyx, 15 CMa and V1449 Aql, the periods identified 
by H2011a and 112011b are in fact those yielding the low- 
est reduced (although in the case of a Pyx there are 
hundreds of comparable choices). In the case of 33 Eri, 
the period selected does not provide the best fit; a longer 
period of 10.437 days yields an improved value of /v. 
Ultimately we see little objective justification for select- 
ing one amongst the various possible period solutions for 



these data, as there are numerous solutions providing an 
acceptable fit to the time series. 

Because the temporal sampling of these data is diverse, 
the data window function is expected to be rather com- 
plex. Therefore to further explore possible periodicities 
in the data, we have performed an independent EFT pe- 
riodogram analysis using the clean procedure as imple- 
mented in CLEAN-NG (Gutierrez-Soto et al., 2009). The 
CLEANed periodograms typically exhibit one dominant 
period, however in only one case (15 CMa) is that domi- 
nant period close to that selected by 112011a and 112011b. 
Moreover, the different reductions of the FORS data do 
not yield the same periods (see Fig. [S]). Application 
of our period analysis procedure to synthetic time se- 
ries with the same window functions and error bars but 
consisting only of noise show that period solutions with 
similar amplitudes to those obtained from the published 
data can be easily obtained for all stars but 15 CMa. 
However, we find that when the most significant pub- 
lished measurements of 15 CMa are removed from the 
analysis (i.e. those with significance larger than 2<t), the 
resulting periodogram continues to show a strong peak 
at that same period, a persistence which can only suggest 
that the reported period arises from a fit to the noise. 

6. CONCLUSIONS 

In this paper we have re-examined the conclusions 
reached by H2011a and H2011b regarding the existence 
of magnetic fields in 6 pulsating B-type stars, using both 
the measurements presented in the original papers and 
re-reduced archival data, as well as our own Stokes V ob- 
servations of these stars obtained using high-resolution 
ESPaDOnS and Narval spectropolarimetry. 

In one case - e Lupi - Least-Squares Deconvolution ap- 
plied to the high-resolution Stokes V spectra confirms the 
presence of a magnetic field, e Lupi therefore represents 
a newly-confirmed magnetic star for which we measure a 
longitudinal field value in very close agreement with the 
mean given by H2011a. 

For the remaining five stars in the sample studied in 
this paper - 15 CMa, 33 Eri, a Pyx, HY Vel and V1449 
Aql - the high-resolution data do not confirm the pres- 
ence of a magnetic field, although in the case of HY Vel 
the precision is quite poor; further observations of this 
star are required to reach any firm conclusions about 
its magnetic properties. For the four stars for which 
magnetic models were reported by H2011a and H2011b 
but for which no signal is detected in the multiple ES- 
PaDOnS and Narval observations - 15 CMa, 33 Eri, a 
Pyx and V1449 Aql - the Stokes V profiles predicted by 
the models are in clear disagreement with our observa- 
tions. An examination of periodograms shows no strong 
reason to prefer the rotational periods derived by H2011a 
and H2011b over a range of other, equally probable pe- 
riods, except perhaps in the case of 15 CMa, although 
we emphasize that for this star there is no indication 
of the predicted structure within the Stokes V profile, 
even for those data which - examining merely the phased 
{Bz) measurements - are in apparent agreement with the 
FORSl/2 data. Furthermore, the periods suggested by 
H2011a are not robust under re-reduction of the FORS 
data. Based on our examination of the published results, 
our re-analysis of the original archival spectra, and anal- 
ysis of new high-resolution spectropolarimetry, we are un- 
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able to confirm the magnetic field detections, periods and 
models presented by H2011a and H2011b for 15 CMa, 33 
Eri, a Pyx and V1449 Aql. 

Many of the inconsistencies between FORSl/2 results 
and those from high-resolution instruments such as ES- 
PaDOnS or Narval can be resolved with the assump- 
tion that the formal error bars of FORS measurements 
have been systematically underestimated by a factor 
of approximately 0.5; alternatively, an increase of the 
FORSl/2 detection threshold from 3ct to 6cr would elim- 
inate the majority of the unconfirmed detections. This 
conclusion is consistent with those reached by Silvester 
et al. (2009) and Bagnulo et al. (2012). 

H2011a point out that many of the stars discussed in 
this paper are both overabundant in nitrogen (e.g., Morel 
et al. 2008) and exhibit pulsations dominated by a non- 
linear radial mode. They propose that the presence of 
a magnetic field may explain these phenomena. Our in- 
ability to confirm the presence of magnetic fields in these 
stars does not support this proposal. Similarly, the large 
equatorial rotational velocities implied by the assumed 
inclination from the line of sight i and the measured 
vsvcii - which are suggested by H2011a to be mostly 
quite high due to near pole-on inclinations for all stars 
but 15 CMa - should be revisited. We ultimately con- 
clude that no evidence exists in either the results pre- 
sented by H2011a or the present work to support a cor- 
relation between photospheric magnetic fields and chem- 



ical peculiarities or pulsational phenomena of SPB and 
/3 Cep stars: magnetism does not appear to be any more 
common in this class of star than in the overall popu- 
lation of B-type stars (approximately 10% of early-type 
stars are detected as magnetic by the MiMeS survey, as 
summarized by Grunhut et al., 2011). This is in contrast 
to results for rapidly oscillating peculiar A-type (roAp) 
stars, for example, in which magnetic fields are consis- 
tently found to be present and related to the magnetic 
field via the oblique pulsator model (e.g. Kurtz et al. 
1997). 
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